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Tj- ■ Abstract 

O" 

f~^ , We performed 2D and full 3D niagnetohydrodynamic simulations of disk accre- 

"^ I tion to a rotating star with an aligned or misaligned dipole magnetic field. We 

Q^' investigated the rotational equilibrium state and derived from simulations the ra- 

Q . tio between two main frequencies: the spin frequency of the star and the orbital 

H I frequency at the inner radius of the disk. In 3D simulations we observed different 

^ ■ features related to the non-axisymmetry of the magnetospheric flow. These features 

may be responsible for high-frequency quasi-periodic oscillations (QPOs). Variabil- 
ity at much lower frequencies may be connected with restructuring of the magnetic 
/\ ' flux threading the inner regions of the disk. Such variability is specifically strong at 

j^ ■ the propeller stage of evolution. 
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1 Rotational Equilibrium State and Two Main Frequencies 



We first investigated slow viscous disk accretion to a rotating star with a dipole 
magnetic field in axisjTnmetric (2.5D) numerical simulations. Quiescent ini- 
tial conditions were developed so that we were able to observe the accretion 
for a long time [1]. These simulations have shown that many predictions of 
the theories developed in the 1970's (e.g., [2]-[4]) are correct. Namely, we ob- 
served that (1) the accreting matter is stopped by the magnetosphere at the 
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Fig. 1. Fluxes in the rotational equilibrium state. M is the matter flux, Lj is the 
angular momentum flux associated with the field (it varies in time but is zero on an 
average), and L^a is an angular momentum flux associated with matter (it is about 
100 times smaller than Lj). 




Fig. 2. Density distribution (background), magnetic field lines (thin lines) and the 
/3 = 1 line (thick line) where the matter ram pressure equals the magnetic pressure. 
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Fig. 3. Left panel: Radial distribution of the angular velocity in the equatorial plane. 
The dashed line shows the Keplerian velocity. Right panel: Density distribution. 

magnetospheric radius r^ ,where the ram pressure pv'^ + p equals to magnetic 
pressure B'^/Sn, and is then lifted up and accretes to the star through funnel 
flows; (2) the magnetosphere influences the structure of the disk to a distance 
r^ Ri (2 — 4)rm; (3) a star may spin-up, spin-down, or may be in the rotational 
equilibrium state depending on the ratio between magnetospheric radius r^ 
and corotation radius rco = {GM/illY^^, as predicted by Ghosh and Lamb 



[3], [4]. We searched for the rotational equihbrium state and investigated this 
state in detail [5] . Namely, we fixed the parameters of the disk and magnetic 
moment of the star, and changed the angular velocity of the star, Q^,. We 
found angular velocity figg at which the angular momentum fiux to the star 
is zero on an average. Figure 1 shows the variation of fiuxes at the surface 
of the star for the run with Tco = 5: the matter fiux M and the angular mo- 
mentum fiuxes Lm and Lf associated with the matter and the magnetic field 
respectively. The fiux Lf (which always dominates over Lm) varies but is ap- 
proximately zero on an average. We found that in the rotational equilibrium 
state Tco/rm ~ lA — 1.7. Figure 2 shows an example of matter fiow in the 
rotational equilibrium state. We see that matter accretes to the star through 
a funnel fiow, and that some field lines are closed and others are infiated or ra- 
dially stretched by the accreting matter. Part of the magnetosphere, however, 
is always closed or only partially open, and this interaction is important for 
the spin-up/spin-down of the star. The fiuxes in Figure 1 vary with time be- 
cause of infiation and reconnection events in the magnetosphere. We observed 
that in the rotational equilibrium state the rotation of the star is locked at 
some value Qeq ~ kQm, where k ~ (1.4 — 1.6) and Q^n is angular velocity 
of the disk at the magnetospheric radius. The number k is smaller (larger) 
at smaller (larger) magnetic moments of the star fi. These two frequencies, 
z/* = (27r)^^r2* and Um = {27T)^^Qm may represent two fundamental frequen- 
cies in each binary system in the equilibrium state. Figure 3 (left panel) shows 
that for our typical case, shown in Figure 2, this ratio z/^ : z/* Ri 3 : 2 (see 
Figure 3, left panel). This ratio is reminiscent of the typical ratio between the 
frequencies of the kHz QPOs observed from many LMXBs [6], [7]. The sur- 
face of the star may be "marked" by hot spots associated with accretion or 
thermonuclear burning at the surface of the star, while the inner disk radius, 
which is approximately equal to magnetospheric radius r^, niay be marked by 
enhanced density in the disk or with enhanced reconnection at the boundary 
between magnetosphere and the disk. Matter often accumulates and forms a 
larger density peak or, ring, at the radius r^ng ~ r^, thus giving a "mark" 
to this region. However, in some cases, the density peak is located at larger 
distances, as is shown in Figure 3 (right panel). The position of this peak 
varies with the accretion rate. It is closer to r^ during periods of enhanced 
accretion. In case of accretion to a misaligned dipole the ring typically breaks 
into two spiral arms (see next section) which rotate with the frequency of the 
inner regions of the disk. 



2 3D Simulations and High-frequency QPOs 



We performed full three-dimensional simulations of the disk-magnetosphere 
interaction [8], [10], [11], for aligned angular momenta of the star and the disk. 




Fig. 4. Magnetospheric flow is different at different density levels: matter blankets 
the entire niagnetosphere at p = 0.2 (left panel), but forms the funnel streams at 
p = 0.4 (right panel). The Q = 30° case is shown. 

and for different misalignment angles between the magnetic and rotational 
axes. We observed that many features of the interaction are similar to those 
in the axisymmetric case. However, in the 3D case, many new features appear 
which are connected with the non-axisymmetry of the flow and which are 
important for analysis of high-frequency QPOs. 

Simulations have shown that matter typically accretes in two streams. In 
cases of relatively high-temperature of the disk [11] or non-stationary accre- 
tion, more than two streams may form. The shape of the magnetospheric flow 
strongly depends on the density. In the same simulation run we see that the 
largest density matter flows in narrow streams, while the lower density mat- 
ter may blanket the magnetosphere completely (see Figure 4). So variability 
may be related to radiation/obscuration by the magnetospheric streams. The 
frequency of this quasi-periodic variability will be between u r^ u^, and the 
frequency of the inner region of the disk, z/^, which is higher than z/* in the 
rotational equilibrium state (see Section 1), but may be lower in some other 
cases. If two or several funnel streams form, then the frequency will be twice 
or several times higher than that expected in case of one stream. 

Other features which may lead to quasi-periodic oscillations are non-axisymmetric 
features observed in the inner regions of the disk. Matter typically forms 
two-armed spiral waves which may contribute to quasi-periodic variability 
(see example in Figure 5). The frequency associated with these features is 
^spirals ^ 2z/m. Wc should uotc that in case of the dipole fleld there are always 
two features, from opposite sides of the accretion disk, so that the frequency 
is twice as the frequency of the disk. 

Signiflcant variability and quasi-periodic variability may be associated with 
the spots at the surface of the star. The shape of the spots, and also distribu- 
tions of density, velocity, matter and energy fluxes, reflect those in the funnel 
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Fig. 5. Example of magnetospheric flow observed in 3D simulations for Q = 30°. 
Inner regions of the disk form a trailing spiral structure. 
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Fig. 6. Density distribution in the hot spots at the surface of the star for different 
B in the relativistic case (from [12]). 
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Fig. 7. Examples of variability curves with relativistic effects (solid lines) and with- 
out relativistic effects (dashed lines) at different G and inclination angles i (from 

[12]). 

streams: at small misalignment angles O < 30°, the spots have the shape of a 
bow, while at large misalignment angles 6 > 75°, they have the shape of an 
elongated bar [8], [12] (see also Figure 6). It is a common belief that the hot 
spots have fixed positions at the surface of the star. Simulations have shown 



that the positions of the spots vary with time. If the disk is relatively cold, 
or the misalignment angle is large (> 45°), then the spots change their shape 
and position only slightly, so that their frequency slightly varies around stellar 
frequency u r^ u^,. However, if the disk is relatively hot and the misalignment 
angle is small (< 30°), then the funnel streams and spots may rotate faster or 
slower than the star, because one of foot-points is dragged by the disk, which 
may rotate faster or slower than the star. This may lead to significant de- 
parture of the associated frequency Ugpot from the z/*. The expected frequency 
lygpot is in between z/* and I'm- Temporary events of precession of the funnel 
streams (and spots) are expected in periods of non-stationary accretion. 



The light curves from the hot spots are different for different misalignment 
angles G and inclination angles i. To investigate these dependencies, we fixed 
hot spots at the surface of the star at some moment corresponding to quasi- 
equilibrium, and then investigated light curves for different O and i assuming 
that all kinetic energy of the flow is converted to black-body radiation [8]. 
We observed that the amplitude of variability may be quite large and ei- 
ther one or two peaks can be observed depending on G and i. Typically, two 
peaks are observed at larger and i. However, in accretion to neutron stars, 
general relativistic, Doppler, and other effects are important (e.g., [13], [14]). 
Typically, the GR effect dominates unless the star rotates extremely fast. We 
performed 3D simulations taking into account GR effects (approximated by 
the Paczynski-Wiita potential) and other effects for neutron stars with param- 
eters close to those in millisecond pulsars. We observed that the shapes of the 
spots are very similar to those in the non-relativistic case [12] (see Figure 6). 
We also observed that light bending effects strongly decrease the amplitude 
of the variability curves [12] (see Figure 7), agreeing with earlier predictions 
(see, e.g. [15]). In extreme cases, like = 90° and i = 15°, both spots are 
observed at the same time, so that the amplitude of variability becomes very 
small (Figure 7, right panel; see also [15]). 



These two main features associated with (1) spots at the star and (2) spirals 
at the inner radius of the disk may lead to high-frequency quasi-periodic os- 
cillations. Simulations have shown that in case of the dipole magnetic field 
there will be two symmetric features (two spirals and two hot spots) and this 
will double the frequency. We should note that there is another important 
frequency, a beat frequency: Vh = ^m ^ ^* which is taken into account in the 
Beat-Frequency and related models (see [16], [17]). This frequency will cor- 
respond to the lower frequency oscillations. If magnetic field of the star is 
dynamically unimportant then the Sonic-Point model may be applicable for 
analysis of QPOs ([18]). However we did not model cases with such a weak 
magnetic field. 
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Fig. 8. Example of accretion at the weak propeller stage: matter accumulates and 
accretes quasi-periodically to the star, forcing magnetic field lines in the magnetic 
"tower" to reconnect (from [24]). 

3 "Propeller" Stage and Low- frequency QPOs 



Simulations often show variability of fluxes at much longer time-scales (see, 
e.g., figure 1). This variability is connected with restructuring of the external 
magnetic field lines in the corona as a result of the disk-magnetosphere inter- 
action. Such variability appears at different rotation rates of the star, but it 
is the strongest at the "propeller" regime of evolution. 



We studied disk accretion to a star in the propeller regime by axisymmetric 
MHD simulations. In the propeller regime the centrifugal force at the external 
boundary of the rotating magnetosphere is larger than the gravitational force, 
so that accreting matter has a tendency to be flung away (e.g., [19] - [21]). We 
observed that there are two qualitatively different regimes: "weak" propeller 
regime (no outflows) or "strong" propeller regime (with outflows). In both 
cases the star strongly spins-down due to the interaction of the magnetosphere 
with the inner regions of the disk. 

The weak propeller occurs when the transport coefficients in the disk (viscosity 
and diffusivity) are relatively small. Both viscosity and diffusivity were incor- 
porated to the code and were considered in the simplified a— approach [22]. 
The weak propeller realizes at relatively small transport coefficients, ttvis < 0.2 
and adif < 0.2. We observed that magnetic field lines inflate, forming a mag- 
netic tower (see also [23]), that the closed magnetosphere expands in the radial 
direction, and matter accumulates near the magnetosphere and is blocked by 
the inflated magnetic flux from accretion. At some point, the accreting matter 
moves forward, forces magnetic field lines in the magnetosphere to recon- 
nect and accretes to the star. Later, the magnetosphere expands again and 
the process repeats [24] (Figure 8). Such quasi-periodic accretion events have 
lower frequencies u^p at larger magnetic moments fi of the star and larger 
n*. The typical frequency of the oscillations lies in the range u^p ^ O.Olz/* to 
Uwp ~ 0.2z/* depending on parameters of the star and the disk. 




Fig. 9. Example of outflows at the "propeller" stage. Background and dashed lines 
show the angular momentum flux carried by matter. Solid lines are magnetic field 
lines. 

The strong propeller (with outflows) was obtained when the viscosity and 
diffusivity in the disk were larger, a^^s > 0.2 and a^a > 0.2. We observed 
that matter is ejected in the form of conical, wide-angle jets as a result of 
magneto-centrifugal forces [26], [27]. The process is quasi-periodic and we ob- 
served many periods of accretion and outflows. An example of outflows is 
shown in Figure 9. Important factors which led to propeller-driven outflows 
are the larger viscosity and diffusivity in the disk, which helped to transport 
momentum from the magnetosphere to the disk. Also, the enhanced viscosity 
helped the disk matter penetrate deeper into the rapidly rotating magneto- 
sphere of the propeller. Both accretion and outflows occur quasi-periodically, 
which is connected with periods of (1) inward viscous accretion; (2) matter 
diffusion through the outer regions of magnetosphere and angular momen- 
tum transport from the magnetosphere to the disk; (3) accretion to the star 
and outburst to the jet; (4) outward motion of the disk (see also [23], [25]). 
The span of frequencies of quasi-periodic oscillations is similar to that in the 
weak propeller case, Ugp ~ (0.01 — 0.2)z/*. Figure 10 shows an example of 
the quasi-variable matter fluxes with a dominant low oscillation frequency 
Ugp ^ (0.015 — 0.02)z/^,. We also noticed that at even higher viscosity, a^is = 0.6, 
the quasi-periodic oscillations became very ordered (see Figure 11). Their fre- 
quency was observed to drift slowly from Ugp ~ O.lz/* at the beginning of the 
simulation to Ugp ~ 0.2z/* at the end of the simulation after > 1000 rotations 
of the star, because the disk moved closer to the star. Similar oscillations were 
observed in a test run with even higher viscosity [28]. 



We should note that the low-frequency variability is also observed in case of 
slowly rotating stars as a result of the restructuring of the magnetic field lines 
threading the inner regions of the disk (see Figure 1). 
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Fig. 10. Matter fluxes to the star and to the jet when Ovis 
Time is measured in periods of Keplerian rotation at r = 1. 



0.3 and a^if = 0.2. 



500 



Matter flux to the star 




520 



Matter flux to the jet 



540 Time 5«' 



580 



600 



Fig. 11. Matter fluxes at Qvis 
Keplerian rotation at r = 1. 

4 Conclusions 



0.6 and adif = 0.2. Time is measured in periods of 



Analysis of the disk-magnetosphere interaction in axisymmetric and full 3D 
MHD simulations have shown that: (1) In the rotational equilibrium state, 
the ratio between frequencies of the inner disk and of the star is z/^ : z^* ~ 
(1.4 — 1.6) for magnetospheric radii r^ ^ (3 — 6)-R*. This ratio may possibly be 
larger for larger r^; (2) Non-axisymmetric features in the inner regions of the 
disk (typically, two-armed spirals, non-axisymmetry of the central density ring 
may lead to QPO oscillations with v ~ 2//^; (3) At small misalignment angles 
G the funnels may precess, and thus may rotate faster/slower than the star; 
the frequency will vary between v^ and v^n or twice these values (depending 
on the number of spots); precession of funnels is expected during periods of 
non- stationary accretion; (4) In all cases the hot spots slightly change their 
position and shape which may lead to low-amplitude QPOs with frequency 
T-' spots ~ ^*; (5) Low-frequency quasi-periodic oscillations were observed at the 
"propeller" stage, with dominant frequency in the range Vp < (0.2 — O.Ol)//^- 
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